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A discrete-state model of the Fi-ATPase molecular motor is developed which describes not only 
the dependences of the rotation and ATP consumption rates on the chemical concentrations of 
ATP, ADP, and inorganic phosphate, but also on mechanical control parameters such as the friction 
coefficient and the external torque. The dependence on these mechanical parameters is given to 
the discrete-state model by fitting its transition rates to the continuous- angle model of P. Gaspard 
and E. Gerritsma [J. Theor. Biol. 247 (2007) 672-686]. This discrete-state model describes the 
behavior of the Fi motor in the regime of tight coupling between mechanical motion and chemical 
reaction. In this way, kinetic and thermodynamic properties of the Fi motor are obtained such as 
the Michaelis-Menten dependence of the rotation and ATP consumption rates on ATP concentration 
\ and its extension in the presence of ADP and Pi, their dependences on friction and external torque, 

as well as the chemical and mechanical thermodynamic efficiencies. 
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I. INTRODUCTION 

\ F Q Fi-ATPase is a ubiquitous protein producing adenosine triphosphate (ATP) in mitochondria [H, Hj]. In vivo, the 
F Q part of this protein is embedded in the inner membrane of mitochondria and is rotating as a turbine when a proton 
current flows through it, across the membrane. This turbine drives the rotation of a shaft inside the hydrophylic 
. Fi part. This latter is composed of three a- and three /3-subunits spatialy alternated as a hexamer (a/3)3 and 
forming a barrel for the rotation of the shaft made of a 7-subunit 0,[|. Upon rotation, the central 7-shaft induces 
conformational changes in the hexamer, leading to the synthesis of ATP in catalytic sites located in each /3-subunit. 
In their experimental work, Kinosita and coworkers have succeeded to build a nanomotor by separating the Fi part 
and attaching an actin filament or a colloidal bead to its 7-shaft. In vitro, ATP hydrolysis drives the rotation of this 
nanomotor, transforming chemical free energy from ATP into the mechanical motion of the 7-shaft 0,0- This motion 
proceeds in steps of 120°, revealing the three-fold symmetry of Fi-ATPase [!, 0, HI- Furthermore, the experiments 
■ have shown that each of these steps is subdivided into two substeps @, [|| ■ The first substep of about 90° is attributed 
to ATP binding to an empty catalytic /3-subunit of Fi, followed by conformational changes of this subunit and of 
^-jl , the whole hexamer as a consequence of its binding to the other subunits. This deformation of the protein generates 
a torque on the 7shaft, hence its rotation by about 90°. The secondary substep of about 30° follows the hydrolysis 
ON ! of ATP and results into the release of the products adenosine diphosphate (ADP) and inorganic phosphate (Pi), 
O ■ completing the 120° stepjl, ©] . 

In our previous work [Q, we carried out a theoretical study of the stochastic chemomechanics of the Fi-ATPase 
molecular motor on the basis of the experimental observations reported in Ref. Q • The stochasticity of the motion is 
the consequence of the nanometric size of the Fi motor making it sensitive to the thermal and molecular fluctuations 
due to the atomic structure of the protein and surrounding medium [ToL fllj . In our paper [l2| , the rotation angle was 
taken as a continuous random variable and the stochastic process ruled by six coupled Fokker-Planck equations for 
the biased diffusion of the angle and the random jumps between the six chemical states considered in our model. This 
continuous- angle model allows us to reproduce accurately the experimental observations of Ref. 0] and, in particular, 
the random trajectories of the 7-shaft with the steps and also the substeps [Hj]. Taking the angle as a continuous 
variable provides a realistic description of the random motion of the motor. 

However, coarser descriptions are often considered in which the angle (or the position in the case of linear motors) 
performs discrete jumps instead of varying continuously. The discretization is considered not only because it is 
always required in order to simulate the random process of diffusion for the angle [l3| , but also because the angle is 
observed to jump by finite amounts corresponding to the steps and substeps observed in the experiments [7|, [8j. This 
observation suggests an alternative description with discrete states associated with the steps and substeps. As long 
as the steps and substeps can be identified with the chemical states of the rotary motor, this fully discrete description 
can be set up on the basis of chemical kinetics. This description could in principle be deduced by coarse graining 
the continuous-angle description, but is most often taken as a model with parameters to be fitted to experimental 
observations. In such coarse-grained descriptions in terms of discrete states, the stochastic process is no longer a 
process with both diffusion and jumps, but becomes a purely jump process ruled by a so-called master equation 
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[lCtlllj. Such discrete-state descriptions provide simpler models of molecular motors, relying on similar but different 
assumptions compared to the finer descriptions using continuous random variables for angle or position beside the 
discrete chemical states. In this context, the question arises of understanding the relationship between the two levels 
of description of the molecular motor. Since both discrete- and continuous-state models of molecular motors have 
been proposed in the literature, it is becoming important to develop methods to relate both levels of description, 
which is the purpose of the present paper. 

The comparison between the two levels of description is interesting not only for the methodology, but also because it 
can reveal important properties of the molecular motor such as the nature of the coupling between their chemistry and 
mechanics, i.e., the question whether the mechanochemical coupling is tight or loose. As discussed in the following, 
this property can be expressed in terms of the nonequilibrium thermodynamics of the molecular motor. 

The paper is organized as follows. In section|TTl we present in detail the discrete-state description of the Fi molecular 
motor. In section lTITl this description is compared with our previous model with a continuous angle variable. In section 
IIV1 we discuss the properties of the Fi motor in the light of the comparison between the discrete- and continuous-state 
descriptions and, in particular, the question of the mechanochemical coupling and its thermodynamic efficiencies. The 
conclusions are drawn in section Ivl 



II. DISCRETE-STATE DESCRIPTION 



A. Chemistry of the Fi motor 

In the discrete-state description, the discrete values of the angle of the 7-shaft correspond to the chemical states so 
that the mechanical motion of the motor is directly controlled by its chemistry. 
The motor is powered by the hydrolysis of ATP into ADP and Pj: 

ATP ^ ADP + Pi (1) 

This reaction is driven by the difference of chemical potential A/i between the three species (ATP, ADP, and P;): 

A^i = /iatp - A^adp - /ip, ( 2 ) 

where the chemical potential fix is equal to the corresponding Gibbs free energy per molecule defined as 

MX =A*x + fcfiTlnH (3) 

with X = ATP, ADP, or Pi, the absolute temperature T, Boltzmann's constant fee = 1-38 10~ 23 J/K, and the reference 
concentration c = 1 mole per liter at which the chemical potential of species X takes its standard value 

The standard Gibbs free energy of ATP hydrolysis takes the value AG = — A/.t = —30.5 kJ/mol = —7.3 kcal/mol 
= — 50 pN nm at the temperature of 23° C, the external pressure of 1 atm, and pH 7 |l4|. We notice that ATP hydrolysis 
provides a significant amount of free energy of —AG = 12.2 k^T above the thermal energy k^T = 4.1 pN nm. Since 
the chemical potential difference © vanishes at equilibrium, the equilibrium concentrations of ATP, ADP, and P; 
satisfy 



[ATP] 



[ADP] [Pi] 



AG 

= exp ~4.910- 6 M _1 (4) 

k B T 

cq 



showing that ATP tends to hydrolyze into its products. In Eq. (U|) and the following, we no longer write the reference 
concentration c , assuming that the concentrations are counted in mole per liter (M). 

The kinetic scheme of our model is based on the phenomenological observations of 120° rotation of the 7-shaft per 
consumed ATP molecule. In accordance with Ref. [7|, the first substep, the 90° rotation of the 7-shaft, is induced by 
the binding of ATP to an empty catalytic site. The second substep, the 30° rotation of the 7-shaft, is induced by the 
release of ADP and Pj. The process can be summarized by the following chemical scheme 

ATP + [0, 7 (0)] W ^ [ATP* , 7(0 + 90°)] ^ [0, 7 (0 + 120°)] +ADP + P 1 (5) 

-> ^ ' W-l ■- ' W-2 < ' 

state 1 state 2 state 1 



From left to right: In state 1, ATP can bind to an empty (0) /3-catalytic site of Fi with the 7-shaft at angular position 
8. State 1 is thus defined by [0, 7(6*)] • Binding of ATP induces the 90° rotation of the 7 -shaft, which we represent by 
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j(8 + 90°) and fills this catalytic site. ATP* stands for any transition state of ATP between the initial triphosphate 
molecule to the products of hydrolysis ADP and P; before the evacuation of the /3-catalytic site. State 2 is thus 
represented by [ATP* ,7(6* + 90°)]. If Fi proceeds to hydrolysis, the products ADP and P; are released together, which 
induces the secondary 30° rotation and empties a /3-subunit. 

Notice that Eq. ([5]) does not necessarily represent the same /3-subunit during different real catalytic states. Following 
our reference data pj], ADP and Pi are released together, but from a different /3-subunit then the ATP binding one. 
In this case, Eq. represents different /3-subunits; one that binds ATP and one that releases ADP and Pi from 
another ATP molecule previously bound to an other /3-catalytic subunit. Furthermore, if the concentration [ATP] 
is below the nanomolar, Fi can perform uni-site catalysis (!a |. In this other case, Eq. ([5]) would represent a single 
/3-subunit in its two main catalytic states. The relevance of Eq. ([5]) to both cases is a consequence of our approach 
which consists in looking only at the contributions of ATP hydrolysis to the rotation of 7. Very recent experimental 
observations have shown that ADP and P; may be released from different /3-subunits Q . If both product molecules 
are released at different times, a more complete model could be considered by adding a third state to the kinetic 
scheme. Here, the scheme is simplified by lumping together the two states subsequent to hydrolysis in a regime where 
the time delay between the releases of ADP and Pi is short enough, as also considered elsewhere (T3 ] . 

Since Fi is a hexamer composed of three /3-subunits, the reactions ([5]) appear three times for the angles 9 + 120°n 
with n = Q, 1,2, so that the motor has a total of six chemical states. However, by the three-fold symmetry of the 
Fi motor, these six states can be regroup three by three since the process repeats itself similarly every 120° in each 
/3-subunit. We can therefore identify the three states of type 1 together in state 1 and the same for state 2, reducing 
the motor dynamics to a process with only two discrete states. 

According to the mass-action law of chemical kinetics, the reaction rates W p in Eq. depend on the molecular 
concentrations in the solution surrounding the motor as follows 

W +1 = fc +1 [ATP] (6) 
W-i = k-i (7) 

W +2 = fc +2 (8) 

W-2 = fc_ 2 [ADP][Pi] (9) 

where the quantities k p (p — ±1, ±2) are the constants of the forward and backward reactions of binding and unbinding 
of ATP or ADP with Pi and, [ATP], [ADP], [Pj] represent the concentrations of each species. is the constant of 
ATP binding often denoted fc on and the ATP unbinding constant k Si while k + 2 is the constant of ATP catalysis 
denoted k cat . The two-state model includes the possibility of ADP and Pi binding and thus ATP synthesis, with 
the corresponding kinetic constant fc_2- Since the rates W p are measured in s _1 , the constant fe+i has the units of 
M _1 s _1 , and fc+2 the units of s _1 , and fc_2 the units of M _2 s _1 . 

An important aspect of the experiments with beads or actin filaments attached to the 7-shaft is that the behavior of 
the molecular motor also depends on the friction £ of the attached objects movin g in the viscous medium surrounding 
Fi, as well as on the external torque r which is applied in some experiments [8l Il6j. As a consequence, the reaction 
constants depend on both the friction ( and the external torque r. Indeed, a discrete-state description in which 

the motion is directly coupled to the chemical reactions leaves no place to the modeling of mechanical aspects such as 
friction and external torque. On the other, a continuous-angle description based on a stochastic Newtonian equation 
of Langevin type would provide a way to include the mechanical forces of friction and external torque by using the laws 
of Newtonian mechanics. This was achieved in our previous paper [l2j , allowing us to study the dependence on friction 
and external torque. In order to take into account these mechanical aspects in the present discrete-state model, we 
have to incorporate the effect of friction ( and external torque r into the reaction constants k p ((, r) (p = ±1, ±2) and 
to find these dependences by comparison with the continuous-angle model [l2| which includes the mechanical aspects. 
This will be achieved in section HTT1 where the discrete-state model will be completed in this way. We notice that a 
discrete-state model with six chemical states has been studied in Ref. [17j in the particular case where the external 
torque is vanishing and for a fixed value of the friction coefficient. Here, our purpose is to investigate situations 
where the molecular motor is submitted to an external torque to understand how the coupling between chemistry and 
mechanics can be formulated in terms of the reaction rates ©-© of the discrete-state model. 



B. Thermodynamics of the Fi motor 

In order for the description to be consistent with thermodynamics, let us summarize here the conditions that the 
chemical and mechanical properties of the motor must satisfy. 

The internal energy of the motor changes with the rotation angle 8 and the numbers N% of molecules (X=ATP, 



4 



ADP, Pi) entering the catalytic sites according to the Gibbs relation 

dE = TdS + rd8 + ^2fi x dN x (10) 
x 

where r denotes the external torque, fix the chemical potential ([3]), T the temperature, and S the thermodynamic 
entropy. The changes in molecular numbers due to the overall reaction {T]) satisfy 

dN ATP = -dN ADP = -dN Pi (11) 

where dNx is the number of molecules of species X entering the motor. Consequently, the Gibbs relation (fT0|) becomes 

dE = TdS + Td8 + Afi dN ATP (12) 

in terms of the chemical potential ([2]) of the overall reaction. In an open system such as a motor, energy and entropy 
vary a priori because of exchanges with the surrounding medium or possibly due to internal variations: 

dE = d c E + diE (13) 
dS = d c S + diS (14) 

The laws of thermodynamics only rule the internal variations: 

first law: d { E = (15) 
second law: d[S > (16) 

Since energy and entropy are state variables, they recover their starting value after completing a cyclic process in a 
motor whereupon their integration over a cycle vanishes, § dE = and § dS = 0. For the isothermal process of a 
nanomotor in the heat bath given by the surrounding medium, the entropy irreversibility produced over a motor cycle 
is thus given by 

j><kS = - i d c S= ^ j(rd6 + AfidN ATP ) > (17) 

which is deduced by using Eqs. (|12[1 - (I16[) for a cyclic process. Introducing the mean angular velocity of the motor in 
revolution per second 

1 /d6\ , , 

and the mean rate of ATP consumption 



2tt \ dt I 



the entropy production is given by 



in terms of the so-called thermodynamics forces or affinities, 2itt/T and A/i/T, and the corresponding fluxes or 
currents, V and R [l8l[l9l. [20j . All these quantities vanish at the thermodynamic equilibrium. Mechanical or chemical 
energies are provided from the exterior if either the torque r or the reaction free energy A[i are non vanishing. 

In general, molecular motors can be driven out of equilibrium by either the external torque r or the molecular 
concentrations of ATP, ADP and Pi entering A/x, possibly in combination. In this regard, the thermodynamic forces 
2ttt/T and Afi/T are independent control parameters. However, in a regime where the chemistry and the mechanics 
of the molecular motor are tightly coupled, three ATP molecules are consumed per revolution which is expressed by 
the conditions that the ATP consumption rate is three times the velocity, i.e., 

tight coupling: V = - R (21) 
3 
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In this case, the entropy production (pOj) becomes 



1 d' S 

tight coupling: l — = AR > (22) 

Kb ere 



with the chcmomcchanical affinity 



2tt t Aii , . 

A EE + —t- 23) 

mechanics chemistry 

which is here written in dimensionless form. In the regime of tight coupling, the mechanical and chemical thermo- 
dynamic forces are thus no longer independent control parameters but are replaced by the unique chemomechanical 
affinity (|23|) . This affinity vanishes at equilibrium. By using Eq. ([3]), we obtain the equilibrium condition for the 
concentrations : 



[ATP] 



[ADP][Pi] 



exp -4= (AG - — r ) ~ 4.9 lO^M" 1 exp (- — -—) (24) 
k B T\ 3 J V 3 k B Tj v 1 



in the presence of the external torque r. This result shows that the chemical equilibrium is displaced in the presence 
of an external torque. In vivo, such an external torque comes from the F D part of ATPase and the transmembrane 
pH difference in mitochondria. This displacement of equilibrium finds its origin in the mechanochemical coupling 
achieved in the Fi rotary motor but happens only if an external torque is enforced on the 7-shaft. 

In the discrete-state description, the discrete values of the angle of the 7-shaft correspond to the chemical states so 
that the mechanical motion is tightly coupled to the chemistry of Fi, a feature which is characteristic of such discrete 
models. Therefore, the relations ([2lT) . (|22|) . and (|23|) apply to such tight-coupling models. 



C. Master equation description 

The chemistry of Fi is a stochastic process because the arrival of each substrate molecule (ATP or ADP with 
Pi) is a random event in time. Accordingly, the time evolution of the molecular motor is described in terms of the 
probabilities P a {t) to find it in one or the other of the two states a = 1,2 of the kinetic scheme (0. These probabilities 
evolve in time because of the random transitions between both states at each reaction. This time evolution is thus 
ruled by the master equation [l(| [ll[ : 

= E [W p (o»iV(i) - W-M^M*)] (25) 

with a sum over the two reactions p = 1,2 and the two states a' — 1,2 before the transition a'^—xj or after the 

reverse transition a — P -*o' ' . The quantity Wp(<7'|(7) is the transition rate per unit time from the state a' to the state a 
due to the reaction p, which can be identified with the rate of the corresponding reaction between the two chemical 
states a and a'. The master equation conserves the total probability P^{t) = 1 for ah times t. Using the specific 
values of the four reaction rates of Eq. ([5]) resumed in Eqs. ([6])-((9|), we write the time evolution of the states 1 and 2 
with respectively empty and occupied catalytic site as 



dPi 
~dT 
dP 2 
dt 



(W-i + W +2 )P 2 - (W +1 + W- 2 )Pi (26) 
(W+i + W-a)Pi - + W +2 )P 2 (27) 



with the normalisation condition Pi + P 2 = 1 always satisfied. 
The mean consumption rates of the different species are given by 



j£ (Aatp) = W + xPi-W-xP 2 (28) 

-^(iVADp) = ^(A^p.) = W^Px ~ W +2 P 2 (29) 
at dt 



6 



Since the 7-shaft rotates by a substep of 90° during the first reaction p = +1 and by a substep of 30° during the 
second reaction p = +2, the mean angle of the 7-shaft evolves in time according to 

4(0) - l {W+iPi - W^P 2 ) + J (W +2 P 2 - W-2P1) (30) 
dt 2 

in radian per second. 

Stationary state. The master equation admits a time-independent stationary solution such that (d/ dt)Po St ^ — 0. 
The stationary solution for each two states is given by 

p(st) _ k -i + k +2 

1 " k +1 [ATP} + k- 1 + k +2 + k- 2 [ADP}[P i ] y ' 
„W = k +1 [ATP]+k^ 2 [ABP][P i ] 

2 k +1 [ATP] + k- 1 + k + 2 + k- 2 [ADP][P i ] [ ' 

where we used the normalization condition P^ + P^ = 1. We notice that this state is stationary in the statistical 
sense because an individual motor continues to fluctuate in time. The probabilities (|31|) - (j32|) define the statistical 
distribution of its fluctuations between the states 1 and 2 after sampling over a long enough lapse of time and over 
many random realizations. 

In the stationary state, the mean angular velocity of the 7-shaft is thus given by 

V=—^(0\ - 1 fc +i fc +2[ ATp ] - fc-ifc- 2 [ADP][Pi] 

2iTdt [>st 3/c+i[ATP] + /j_i+£; + 2 + fc-2[ADP][Pi] 1 ' 

and the mean rates of molecular consumption by 

R = j t (N ATP ) st = -A(iv AD p) st = -j t (N Pi ) st = W (34) 

The condition (j2~l|) of tight chemomechanical coupling is thus satisfied. 

Equilibrium state. At thermodynamical equilibrium, all the noncquilibrium constraints vanish. In this case, the 
stationary solution represents the equilibrium probability pj cq ^ and obeys the detailed balance conditions [2l[ 

W p {a'\a)P^ =W- p {aW)P^ ] (35) 

for all reactions p = 1,2 and all the transitions a, a' = 1,2. As a consequence, the mean velocity (|33f and the rate 
(J3"4"| vanish which occurs under the condition: 



[ATP] 



[ADP][Pi] 



-r±r± (36) 



Comparing with the condition obtained from the thermodynamics of the motor, we get the equilibrium constant of 
the reaction: 

which depends on the external torque r acting on the molecular motor. We see that considerations of equilibrium 
thermodynamics give a constraint allowing us to fix one reaction constant if we know the three other constants. 
However, equilibrium thermodynamics does not provide enough relations in order to determine the four reaction 
constants. We have therefore to use data from the kinetics of the motor, i.e., when the motor is out of equilibrium. 
This is the purpose of the following subsection. 



D. Determination of the reaction constants 



The motor is in a nonequilibrium stationary state if the chemomechanical affinity (|23[) is non- vanishing, i.e., if the 
equilibrium condition (|36D is not satisfied. In this case, the mean velocity (|33p does not vanish and can be compared 
with experimental data 0] or with the continuous-angle model [l2| which is complete with respect to the mechanics. 
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In this way, each one of the reaction constants k p can be determined by appropriate limits at fixed values of £ and r, 
as explained here below. 

First of all, the mean velocity ()33|1 can be rewritten in the following form: 

v = K, ax ([ATP]-g oq [ADP][Pi]) 

[ATP] + Km + -Kp[ADP][Pi] [M) 



in terms of the constants 



Knax = £ fc+2 (39) 

i^M - (40) 

Kp 55 F 1 (41) 



together with the equilibrium constant i"T cq defined by Eq. (|37|) . The knowledge of these four constants is equivalent 
to knowing the four reaction constants k±i and k±2- Each one of the three constants (|39p. (l4"0"l) , and l|4ip can be 
determined in a specific regime of functioning of the molecular motor. 

The constants ([39| and (|40|) can in particular be determined in the absence of the products of hydrolysis, [ADP][Pi] = 
0, in which case the velocity follows a typical Michaelis-Menten kinetics: 

= WATP] 

V -[ATP] + K M (42j 

Here, we see that Eq . ([40]) is the Michaelis-Menten constant defined as the ATP concentration at which the velocity 
V equals Kn ax /2 [22|. This constant characterizes the crossover between the regime [ATP] <C Km where the rotation 
is limited by the slow arrival of ATP molecules and the saturation regime [ATP] 3> Km where the velocity reaches 
its maximum value fixed by the finite rate of release of ADP and Pi. Accordingly, the constant ()39j) is the maximum 
angular velocity of the motor, which is reached in a regime where the ATP concentration is large enough with respect 
to the Michaelis-Menten constant: 

V max = lim V (43) 

[ATP]->oo 

the limit meaning that [ATP] ^> Km- The Michaelis-Menten constant itself can be determined at low ATP concen- 
trations according to 



this other limit meaning that [ATP] <C Km- 

The third constant (|41| characterizes the decrease of the velocity as the concentrations of the products ADP and P; 
are increased. It turns out as we shall see in the following that this constant is larger by several orders of magnitude 
than the equilibrium constant: 

K P » K cq (45) 

Accordingly, the term involving the equilibrium constant in the numerator of Eq. (|38[) can be neglected. Hence, the 
constant Kp can be determined in the regime where 

[ ATP ] + ^«[ADP k2 [ Pi ]«ra (46) 



K P J K c ^ 

by taking the difference of the inverses of the velocities at two different concentrations of ADP, keeping fixed the other 
concentrations : 

[ATP] A(l/V) 

K P -V max [p;] A[ADp] (47) 



with the notation AX = — X\ . 
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Another consequence of the inequality (|45[> is obtained by using the definition (|37p of the equilibrium constant: 

fc -i = IF* fc +2 « ^+2 (48) 
whereupon the Michaelis-Menten constant (|4"0)l is essentially independent of in the present case: 

Km — -p— (49) 

K+l 

and thus directly determines the constant k+\ of ATP binding once the constant k+2 of product release is obtained 
thanks to the maximum velocity (I39|) . Subsequently, the constant k-2 of the binding of ADP and Pi is determined 
from the constant Kp as explained here above. Finally, the constant k-i of ATP unbinding is determined with the 
equilibrium constant ([57)1 . In summary, we find successively: 

(50) 
(51) 
(52) 

K cq (53) 



k+2 


— 3 V^nax 


k+\ 


k+2 
K M 


k-2 


= k+iKp 


k-i 


k+ik+2 

k-2 



III. COMPARISON WITH THE CONTINUOUS-STATE DESCRIPTION 



In this section, we compare the aforementioned discrete-state description with a continuous-state description we 
previously reported on |l2j. The continuous-state description considers the rotation angle as a continuous random 
variable instead of supposing the angle to jump by a finite amount at each reactive event. Accordingly, the continuous- 
state description allows us to incorporate mechanical aspects in the modeling which should otherwise be assumed in 
a discrete-state model, as explained here above. Therefore, the continuous-state description provides the dependence 
of the rotational velocity on both the friction £ and the external torque r and is thus more complete in this regard. 
Consequently, the dependence of the reaction constants of the discrete-state model on friction and external torque 
can be determined by comparison with the results of the simulation of our continuous-state model [l2j , which is the 
purpose of the present section. 



A. Fokker-Planck equation description 

In the continuous-state model the system is found at a given time t in one of the six chemical states a = 1, 2, 6 
and the 7-shaft at an angle < 8 < 2ir. There are six chemical states because the three /3-subunits can be either 
empty or occupied by a molecule of ATP or by the products ADP and Pi of hydrolysis. 

Consequently, the system is described by six probability densities p a {9 1 t) normalized according to 
Y^j=i Jo ^ P<y(9, t)d6 — 1. The time evolution of the probability densities is ruled by a set of six Fokker-Planck 
equations coupled together by the terms describing the random jumps between the chemical states a due to the two 
chemical reactions (ATP binding and release of the products ADP and P^) with their corresponding reversed reactions 

d tP(T (9,t) + de J<r(6,t)= 53 J2 bv(M)v^(«)-P«r(M)M- ft ^v(9)] (54) 

p=l,2< T '(^<j) 

where the probability current densities are given by 

J„(6,t) = —D d eP a(9,t) + i [-doU a {6)+T]pa(6,t) (55) 

The diffusion coefficient D can be expressed in terms of the friction coefficient £ according to Einstein's relation 
D = kftT/C,. The friction coefficient £ can be evaluated for a bead attached to the 7-shaft, a bead duplex, or a 
cylindrical filament, as presented elsewhere @, HH, [fl HH] ■ In the case of a bead of radius r attached off-axis with its 
distance at a distance x — r sin a from the rotation axis, the friction coefficient is given by 



C = 27r?7r 3 (4 + 3 sin 2 a) 



(56) 
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with the water viscosity rj = 10~ 9 pN s nm -2 and a = ir/6 fl2l [T3|. 

When the motor is in the chemical state a, the 7-shaft is submitted to the external torque r and the internal torque 
—dgU a due to the free-energy potential U a {9) of the motor with its 7-shaft at the angle 9. Applying an external 
torque to the motor has the effect of tilting the potentials into U a {9) — t9 which eases the rotation or makes it harder, 
depending on the sign of r. 

B. Coarse graining into a discrete-state model 

The correspondence with the discrete-state description in terms of the master equation (|25|) can in principle be 
established by coarse graining the continuous angle into discrete states. These discrete states correspond to the angular 
intervals 9 a < 9 < 9 a + 2tt/3 where the 7-shaft spends most of its time while in the chemical state a. Accordingly, 
the probabilities ruled by the master equation (|25[) are related to the probability densities of the continuous-state 
description (|54|) according to 

/■e <r +27r/3 

P*(t)= / Pa(9,t)d9 (57) 

where the angular integral is carried out over the aforementioned intervals. In this way, a fully discrete description 
could be inferred from the continuous-state description. In general, this reduction from one description to the other by 
the aforementioned coarse graining leads to non-Markovian equations. In the case where there is a net separation of 
time scales between the dwell times and the transit times between the discrete states, the non-Markovian effects may 
be negligible and a description in terms of a Markovian equation such as the master equation P5"|) may be obtained. 
This is the situation we here consider. 

In this case, the transition rates W p (a'\cr) of the master equation ([25| can be deduced from the solution of the 
Fokker-Planck equations (f54"|) . We emphasize that the transition rates W^ct'Ict) of the master equation ([25]) do not 
take the same values as the rates w p _ a '- iCr (9) appearing in the Fokker-Planck equations (|B"4")) . Indeed, the latter ones 
depend on the angle although the former ones do not. A method of deduction of the former from the latter with Eq. 
(|57p would show that the transitions rates of the discrete-state description are given by intergrals of the continuous- 
angle rates in Eq. (|54p combined with the probability densities of the quasi-stationary dwell states corresponding to 
the substeps. In spite of their conceptual and numerical differences, the transition rates appearing in both descriptions 
are in correspondence as far as the chemical reaction they describe are concerned, as shown in Table [U 

TABLE I: Comparison of the transition rates of the discrete model ruled by the master equation (I25|l with those of the continuous 
model ruled by the Fokker-Planck equations (|54[1 . The transition rates are in correspondence by the chemical reaction they 
describe: ATP binding for p — +1; ATP unbinding for p — — 1; the release of ADP and Pi for p — +2; the binding of ADP and 
Pi for p = —2. We notice that the rates of the continuous model depend on the angle 9 of the 7-shaft. This dependence is of 
Arrhenius type in terms of the inverse temperature f3 = (feT) -1 and the free-energy potentials U(8) and U{9) (for respectively 
the empty and occupied catalytic sites) and [/* (9) and (7* (9) (for respectively the transition states of ATP binding or unbinding 
and of the release or binding of ADP and Pi). See Ref. [131 f° r the forms of these potentials and the numerical values of the 
parameters of our continuous model. 



W p (a'\a) 


W p ^l^ t y(9) 


fc+i[ATP] 

k-i 

k+2 
fc_ 2 [ADP][Pi] 


fc [ATP] exp {-/? [17* (0) - U{9) - G° ATP ] } 
fc exp{-/3 [[/*(£)- (7(0)]} 

fcoexp{-/3[t/*(0)-t/(0+2i)]} 

feo[ADP][Pi]exp{-/3|[7t(0) -U(9)- G° ADP - G° P ]} 



We notice that the rates w p . a i appearing in the Fokker-Planck equations (|54f only concern the chemical 
reactions and, therefore, do not depend on the friction coefficient Q and the external torque r which only enter in the 
current densities (|55[) appearing in the left-hand side of the Fokker-Planck equations (|54[) . In contrast, the reaction 
constants k p shown in Tablc[T]nccessarily depend on the friction coefficient £ and the external torque r, which otherwise 
would not appear in the discrete model given by the master equation (|25[) . Accordingly, we now need to determine 
the dependences k p {Q,T) of the discrete-state reaction constants on both friction and external torque, which is the 
purpose of the next subsection. 
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TABLE II: Values of the coefficients of the Taylor expansions (|59[) - (|60[ ) of the functions a p (r) and b p (r) giving the reaction 
constants fc+i(C, 7"), k+i{C,,T), and fe_2(f, r), according to Eq. (158[l . 



coefficient 


fc+l(C,T) 


fc+2(C,r) 


*-a(C,T) 


units 






-16.952 


-5.973 


-19.382 






«« 


9.8 10~ 4 


1.7 10~ 4 


1.29 10 _1 


(pN nm) 


-l 


o, W 


5.8 1(T 4 


1.0 10" 3 


2.8 10~ 4 


(pN nm) 


-2 


6 (o) 

Op 


-16.352 


-2.960 


-18.338 






Op 


-6.6 10" 2 


-2.7 10~ 2 


5.9 10" 3 


(pN nm) 


-1 


b (2) 

Op 


1.0 10" 3 


3.6 10" 4 


-2.1 10" 4 


(pN nm) 


-2 



C. Dependence of the reaction constants on friction and external torque 

The original and key feature of our work is the inclusion in the discrete-model reaction constants k p (p = ±1, ±2) 
of the mathematical dependence on friction and external torque. This dependence expresses the coupling between the 
chemical and mechanical properties of the motor. 

We determine the dependences of the reaction constants k p on both friction and external troque by fitting them to 
the simulations of our continuous-angle model [13] • Since this latter has been fitted to experimental data, the present 
fitting procedure is comparable to a fitting to the experimental data of Ref. Q. We use the method presented in 
subsection III Dl In this way, we evaluate successively the reaction constants thanks to Eqs. (|50 p - (|53|) for different 
values of the external torque r and the friction £. 

In this regard, an essential observation is that the reaction constants become independent of the friction coefficient 
C at low friction and decrease as the inverse of the friction at high friction: k p oc C, . The motor is functioning 
in a reaction-limited regime at low friction and in a friction-limited regime at high friction (T^ |. At high friction, 
the substeps are no longer visible since the rotation is slow down by the friction, in which case the Fokker-Planck 
equations (|54|) suggest indeed that the rate constants should scale as 

The crossover between the low- and high-friction regimes can be well described by giving the following analytical 
form to the reaction constants: 

MC,r) = eBp(T)+ 1 etp(T)c (58) 

with p = ±1, ±2. The coefficients of the function in the denominator are taken as exponentials in order to guarantee 
the positivity of the reaction constants, the friction coefficient ( being always non- negative. The coefficients e ap ' T ) 
can be determined in the low- friction regime and the coeffcients e bp ^ in the high- friction regime. 

The dependence on the external torque r appears in the functions cl p {t) and 6 p (t), which are taken as expansions 
in powers of r limited to the second order: 

a P (r) = aW+a^T + a^S + OiT 3 ) (59) 
bp {r) = b^+b^r + b^r 2 + 0(r 3 ) (60) 

with p = 1,2. The coefficients of these expansions are fitted in intervals of values of the external torque which 
are typically |t| < 20 pN nm. The values of the coefficients of Eqs. (|59"|) - (f6"0"|) for the reaction constants k + i((,r), 
A;+2(C) t )j an d fc-2(Ci r ) are given in Table ITu 

The last constant for ATP unbinding is finally obtained by using Eq. (|53|) as 

with AG = -50 pN nm. 



IV. PROPERTIES OF THE Fi MOTOR 



The comparison between the discrete and continuous descriptions sheds a new light on the properties of the Fx 
motor. On the one hand, the continuous- angle model reproduces the experimental observations of Ref. Q and its 
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simulation can test the assumptions of the discrete-state model, in particular, the assumption of tight coupling. On 
the other hand, the discrete-state model can be treated analytically thanks to its simplicity. In this way, the fitting 
of the discrete-state model to the continuous one allows us to investigate more closely the properties of the Fi motor 
in the regime of validity of the discrete-state model. 

A. Tight versus loose chemomechanical coupling 

In order to determine the regime of tight coupling between the chemistry and the mechanics of the Fi motor, both 
the angular velocity V and the ATP consumption rate R have been simulated with the continuous-angle model (|54| for 
different values of the external torque r and chemical potential difference A/x, which are the corresponding affinities. 
We depict in Fig. Q]the curves in the plane (r, A/i) where either the velocity vanishes, 1^ = 0, or the ATP consumption 
rate vanishes, R = 0. The value of the external torque where the velocity vanishes is the so-called stalling torque. 
Since the variations of the chemical potential difference A/i can only be one dimensional in such a representation, we 
take a combined variation of the concentrations of the difference species as explained in the caption of Fig. [TJ The 
two curves V = and R = intersect at the origin (r = 0, A/i = 0) which is the thermodynamic equilibrium point. 
We notice that the curve V = is above the curve R = in the plane of the chemical potential difference A/z versus 
the torque as it should according to the second law of thermodynamics ([20]) . 




torque (pN nm) 

FIG. 1: Chemical potential difference A/i in units of fcaTln 10 versus the external torque r for the situations where the rotation 
rate V (circles) and the ATP consumption rate R (squares) vanish in the continuous model (|54| and compared with the 
straight line A/i = —2nr/3 where the chemomechanical affinity ([23} vanishes, A = 0. The concentrations are fixed according 
to [ATP] = 4.9 lO 8 "- 11 M and [ADP][Pi] = io" - 20 - 5 M, in terms of the quantity a = A/x/(fc B Tln 10). The bead attached 
to the 7-shaft has the diameter d — 2r — 80 nm and the temperature is of 23 degrees Celsius. The torque where V — is 
called the stall torque. The determination of the curves V = and R — is difficult close to the thermodynamic equilibrium 
point (t = 0, A/i = 0) because both the rotation rate V and the ATP consumption rate R are very small in this region, which 
explains the absence of dots close to the origin. 

In the tight-coupling regime, the condition (|2ip should hold, which implies that the mechanical and chemical 
affinities r and A/i are no longer independent but should combined into the unique chemomechanical affinity (|23[) . In 
this regime, the vanishing of the rates, V — and R — 0, should thus occur on the curve where the chemomechanical 
affinity (|23p vanishes, i.e., along the straight line 

A^ - -— r (62) 

This is observed in Fig. [T] for values of the external torque extending from zero down to about r ~ — 30 pN nm 
and for chemical potential difference from zero up to Afi ~ 14 k^T, which delimits the zone where the tight-coupling 
assumption is satisfied. Outside this zone for higher values of A/i or lower values of r, the Fi motor is no longer 
functioning in the tight-coupling regime and the chemomechanical coupling becomes loose [24| . 
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Since the coincidence of the curves V — and R — along the straight line (|62|) is a feature of the discrete-state 
model (|27p . we may expect it describes the Fx motor in the aforementioned zone of tight coupling. 



B. Rotation rate versus ATP concentration 

Random trajectories of the discrete model can be simulated thanks to Gillespie's numerical algorithm (2f| [2(1 . 
Examples of random trajectories are depicted in Fig. [2] for different values of ATP concentration, illustrating the 
Michaclis-Mcnten kinetics described by Eq. (|42[) . At low concentrations [ATP] <C Km — 17/JVI, the motor is 
essentially waiting for the arrival of new ATP molecules with its 7-shaft at some angle 120°n (n integer). Instead, 
at high concentrations [ATP] 3> Km, the rotation is limited by the release of the hydrolysis products ADP and P; 
from its catalytic sites with its 7-shaft at some angle 90° + 120°n (n integer), as clearly seen in Fig. [2j In this 
respect, the random trajectories of the discrete model reproduce the jumps by 90° and 30° corresponding to the 
substeps experimentally observed in Ref. 0]. We notice, however, that the discrete model cannot reproduce the 
small- amplitude fluctuations around the dwell angles as the continuous-angle model does [l2l |. 

The crossover between the regime at low ATP concentration and the saturation regime of the Michaelis-Menten 
kinetics is seen in Fig. [3] where we directly compare Eq. (|4"2")l of the discrete model with experimental data from Ref. 
Q. At low ATP concentrations, the rotation rate is proportional to the ATP concentration while the rotation rate 
reaches its maximum value of about 130 rev/s in the saturation regime. 

In order to appreciate the nonequilibrium thermodynamics of the molecular motor, it is interesting to depict the 
rotation rate as a function of the affinity (|23|) instead of the ATP concentration. Indeed, the former is a substitute 
of the latter, as shown by expressing the concentrations in terms of the chemical potentials by Eq. ([3]) and using the 
definitions of the chemomechanical affinity (|23| and of the equilibrium constant (|3T[) to get 

[ATP] = K cq [ADP] [Pi] e A (63) 

In this way, we recover the equilibrium relation (|36[) between the concentrations with Eq. (|37p in the thermodynamic 
equilibrium state A = corresponding to given concentrations of ADP and Pj. Substituting Eq. (|6"3"]) into Eq. (|38p . 
we obtain the following expression for the rotation rate: 

V = V r\ (CA sv2 ( 64 ) 



with the constant 



= 3IW^ cq [ADP][Pj] ( 
~ (tfeq + ATpHADPlpil+tfM 1 j 



This coefficient controls the linear response of the molecular motor because 



V* \\LA^A«l 
I/ max for A > 1 



The analytic form ([B^]) shows that the rotation rate depends on the thermodynamic force A in a highly nonlinear way, 
in contrast to what is often supposed. The nonlinear dependence is very important as observed in Fig. [4] The linear 
regime extends around the thermodynamic equilibrium point at Afi = where the function V(A) is essentially flat 
because the linear-response coefficient takes the very small value L ~ 10 s -1 . Since the affinity is about A ~ 21.4 
under the physiological conditions [ATP] ~ 1(T 3 M, [ADP] ~ 1CT 4 M, and [Pi] ~ 1(T 3 M *14], the rotation rate 
would take the extremely low value V ~ LA/3 ~ 6.5 rev/day if the motor was functioning in the linear regime. 
Remarkably, the nonlinear dependence of Eq. (164|) on the affinity A allows the rotation rate to reach the maximum 
value Vi n ax — 130 rev/s under physiological conditions. 



C. Rotation rate versus friction 



In Fig. [5l we show the effect of friction on the motor's velocity in the absence of ADP or Pi. At low friction, the 
velocity saturates exhibing the reaction-limited regime. At high friction, we see the rapid decrease of the velocity 
with increasing friction in the friction-limited regime [l2j]. In this latter regime, the velocity decreases as the inverse 
of the friction coefficient V cx C _1 in consistency with the analytical form (|5"8|) given to the reaction constants. A 
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0.02 0.04 0.06 0.08 0.1 

time (s) 

FIG. 2: Simulation of random trajectories of the discrete-state model for [ADP][Pi]=0, a temperature of 23 degrees Celsius, a 
bead of diameter d — 2r = 40 nm, and a zero external torque. 



1000 




[ATP] (M) 



FIG. 3: Mean rotation rate of the 7-shaft of Fi in revolutions per second, versus the ATP concentration [ATP] in mole per liter 
for [ADP][Pj]=0. In accordance with the experimental setup pj, the diameter of the bead is d = 2r = 40 nm, the temperature 
is of 23 degrees Celsius, and the external torque is zero. The circles are the experimental data from Ref. 0. The solid line is 
the result of numerical simulation of the two-state model. 




FIG. 4: Mean rotation rate versus the chemical potential difference in units of the thermal energy faT. The thermodynamic 
equilibrium corresponds to Afi = 0. The ATP concentration is given in terms of the chemical potential difference by [ATP] = 
[ADP][Pi]exp[(A^i - A/i )/(teT)] ~ 4.9 10~ 6 M" 1 [ADP] [Pi] exp[A/x/ (fcsT)] since A/j. = -AG = 50 pN nm. The results of 
the discrete model (solid lines) are compared with the continuous model (dots) for three different values of [ADP][Pj]. The 
diameter of the bead is d = 2r = 40 nm, the temperature 23 degrees Celsius, and the external torque zero. 
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1000 




viscous friction (pN nm s) 



FIG. 5: The mean rotation rate in revolutions per second versus the viscous friction coeffcient £ for a bead, a bead duplex, or 
a filament attached to the 7-shaft. The circles and triangles are the experimental data in Fig. 2 of Ref. at [ATP] = 2 mM 
(squares and downward triangles) and [ATP] = 2 fj,M (circles and upward triangles). The squares and circles correspond to the 
single beads and bead duplexes, the triangles to the actin filaments. The solid lines are the results of the present model with 
[ADP][Pi] = 0, the temperature of 23 degrees Celsius, and zero external torque. 




FIG. 6: Mean rotation rate of a bead of radius r — 20 nm attached to the 7-shaft in revolutions per second versus [ADP] 
(in mole per liter) and varying ATP concentrations. The concentration of Pi is 10 -3 M and the temperature of 23 degrees 
Celsius, and the external torque vanishing. The simulations of the continuous- angle model (squares, diamonds, and circles) are 
compared with the analytical expression (|38p for the discrete model of Table HT1 (lines). 



similar behavior is observed at [ATP] = 2 mM in the saturation regime of the Michaelis-Menten kinetics and at ATP 
concentrations lower than the Michaelis-Menten constant for [ATP] = 2 /iM < Km = 17 /iM. The agreement is as 
good as for the continuous- angle model 



D. Rotation in the presence of ATP hydrolysis products 



In the presence of ADP and P, in the environment of the motor, Eq. (|38p shows that the rotation rate decreases, 
as expected since these products tend to counteract ATP hydrolysis that is powering the motor. This phenomenon is 
known as ADP inhibition 0, |H| • There are two possible causes of the decrease of rotation rate if the concentrations 
of ADP and Pj are positive: (1) the term in the numerator of Eq. (|38|) where the equilibrium constant K eq multiplies 
[ADP] [P;] ; (2) the term in the denominator of (f38|) where the constant Kp multiplies [ADP] [Pi]. As explained in 
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(b) [ATP] =4.9 10" 6 M 
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FIG. 7: Mean rotation rate V of the 7-shaft of Fi in revolutions per second versus the external torque for (a) [ATP] == 4.9 10 7 M 
and (b) [ATP] = 4.9 10 -6 M (circles for the continuous-state model and solid line for the discrete-state model). The other 
concentrations are fixed to [ADP] = 10~ 4 M and [Pi] = 1CT 3 M. The squares show the reaction rate R divided by three in order 
to display the regime of tight coupling where V = R/3. The diameter of the bead is d = 2r = 160 nm and the temperature of 
23 degrees Celsius. 



subsection III Dl and can be checked with Table HT1 the inequality (j¥5|) holds by several orders of magnitude so that the 
main cause of the decrease of the rotation rate is the term in the denominator of Eq. (|38p due to the reaction p = — 2 
of binding of ADP and P; to the catalytic sites of the Fi motor. 

As observed in Fig. [6]which compares the continuous and discrete models, this effect manifests itself above millimolar 
concentrations of ADP if inorganic phosphate is in millimolar concentration. We notice that the decrease of the 
rotation rate goes as the inverse of the concentrations of ADP and Pi, V cx ([ADP] [Pi]) -1 as described by Eq. (f35|) 
given the fact that the inequality (I45|) holds. 



E. Dependence on the external torque 

Figure [7] depicts the dependence of the rotation and ATP consumption rates on the external torque r for both the 
continuous-angle and discrete-state models, showing their agreement in the range of validity of tight coupling. This 
range of validity is observed in Fig. [T]to extend down to r ~ —30 pN nm. The reason is that the mechanical motion 
becomes decoupled from the chemical reactions if the external torque is too much tilting the internal free-energy 
potential surfaces of the motor. For vanishing or moderate values of the external torque, the potential surfaces have 
barriers which stop mechanical motion and allow for jumps between potential surfaces corresponding to different 
chemical states of the molecular motor, coupling in this way the mechanical motion to the chemical reactions. Beyond 
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some threshold for the external torque, these barriers disappear causing the decoupling between mechanics and 
chemistry (l2j . This decoupling may happen at both negative and positive thresholds for the external torque, as seen 
in Fig. [7] for the continuous-angle model. The tight-coupling condition V = i?/3 is only observed in the interval 
|t| < 30 pN nm for the continuous-angle model, although it always holds by assumption for the discrete-state model. 
Therefore, the range of validity of the discrete model is restricted to the interval |r| < 30 pN nm of values of the 
external torque, as already discussed about Fig. [T] 

Figure [7J also shows the remarkable feature of the discrete model to reproduce the phenomenon of stalling torque 
that the mean rotation rate can be stopped at a negative critical value of the external torque opposing the rotational 
motion. Indeed, the value of the external torque where both the rotation and ATP consumption rates vanish can 
be obtained in the discrete model by Eq. ([M]) . which shows that the chemomechanical affinity should vanish at the 
stalling torque. Whereupon, we recover the condition (|6l?|) giving the stalling torque as 

T = -i_A// = -—( A M ° + fc B rin J ATP] \ (67) 
sta11 2tt p 2tt V [ADP][Pi]y V ' 

in the tight-coupling regime. For [ADP][P;] = 10 _7 M 2 , the chemical potential difference takes the values A/i = 
56.5 pN nm and A/i = 65.9 pN nm, for respectively [ATP] = 4.9 10" 7 M and [ATP] = 4.9 10~ 6 M. Hence, Eq. (j67|) 
gives respectively the stalling torques r s tali = —27.0 pN nm and r s t a ii = —31.5 pN nm, as indeed observed in Figs. [7£i 
and [7b. 



F. Chemical and mechanical efficiencies 



Under a negative external torque r < 0, the Fi motor can synthesize ATP, in which case the ATP consumption rate 
as well as the rotation rate are negative, R < and V < 0. In this regime of ATP synthesis, the chemical efficiency is 
defined as the ratio of the free energy stored in the synthesized ATP over the mechanical power due to the external 
torque: 

Vc = -7. 77 (68) 

27TT V 

such that < ?7 C < 1 (27[. 

A mechanical efficiency can similarly be defined in the regime where the rotation is powered by ATP as the inverse 
of the chemical efficiency [27j : 

*» S -A£fl (69) 

The mechanical efficiency satisfies < r] m < 1 in the regime where the external torque is still non-positive while both 
the rotation rate and the ATP consumption rates are positive, V > and R > 0. 

It is known [27| that the chemical and mechanical efficiencies reach their maximum values under different conditions 
as shown in the Appendix. In the tight-coupling regime where R = 3V, these conditions coincide and the chemical 
and mechanical efficiencies ([6"5]l and (|6"9")l become 

. , n . 1 3Au 
tight coupling: r/ c = — = (70) 

T)m 27TT 

in agreement with Eq. (|A27[) derived in the Appendix from the assumption of linear response. In the tight-coupling 
regime, the chemical and mechanical efficiencies can reach the maximal unit value at the stalling torque where Eq. 
([6"Tj) holds. This remarkable result is observed in Fig. [5] depicting the chemical and mechanical efficiencies versus the 
external torque under conditions corresponding to the chemical potential difference A/i = 56.5 pN nm. 

In the interval r sta ii = —27.0 pN nm < r < 0, the external torque is opposed to the mean rotation rate but the 
motor consumes ATP, so that the mechanical efficiency (|69l) is postive. In this interval, the coupling is tight so that 
the mechanical efficiency computed with the continuous-angle model (squares) agrees with the prediction (|70p of tight 
coupling which holds for the discrete model. Indeed, the mechanical efficiency reaches the unit value at the stalling 
torque marked by the vertical line. 

On the other hand, the motor synthesizes ATP under the action of the external torque for r < T sta \\ = —27.0 pN nm, 
where the velocity becomes negative and the chemical efficiency (f6"5)l positive. As explained above, the motor is no 
longer in a regime of tight coupling below the stalling torque so that the chemical efficiency computed with the 
continuous-angle model differs from the prediction ([701) °f tight coupling. It is only at the stalling torque that both 
agree while they reach the unit value. However, as seen in Fig. [51 the actual chemical efficiency is lower than possible 
by tight coupling because of the loose chemomechanical coupling below the stalling torque. 
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FIG. 8: Chemical efficiency (I68|l and mechanical efficiency (I69|l versus the external torque r in the continuous-state model 
(respectively circles and squares joined by a solid line) and compared with the prediction (|70[) of tight coupling (dashed 
lines). The vertical solid line indicates the stalling torque at r = —27.0 pN nm. The concentrations are [ATP] = 4.9 10~ 7 M, 
[ADP] = 10~ 4 M, and [P ; ] = 10" 3 M. The diameter of the bead is d = 2r = 160 nm and the temperature of 23 degrees Celsius. 
These conditions are identical as in Fig. [7^. The predictions of tight coupling (dashed lines) are respectively rj c — r sta n/r for 
T < T sta n, and 77 m = r/r sta n for T sta ii < r < 0, with r sta n = -27.0 pN nm. 



V. CONCLUSIONS 



In the present paper, we have studied one of the simplest possible stochastic processes describing the stochastic 
chemomechanics of the Fi-ATPase molecular motor, as observed in Ref. This description considers the two 

discrete states corresponding to the steps and substeps observed in the rotary motion of the 7-shaft of Fi • This 
two-state discrete model is based on the master equation ruling the time evolution of the probabilities of the two 
discrete states. The model is analytically solvable, which provides us with an interesting insight in our understanding 
of the behavior of the Fi motor. The discrete model is set up by fitting its rate coefficients to simulations with our 
previously continuous-angle model [l2[ ■ Since this latter was itself fitted to the observations of Ref. 0] , both the 
discrete and continuous models describe the same experimental system. 

The comparison between the discrete-state and continuous-angle descriptions reveals important properties of the Fi 
motor. Indeed, the discrete-state model presupposes that the mechanical motion and the chemical reactions powering 
the motor are tightly coupled, although the continuous-angle model does not. Accordingly, the comparison between 
both models clearly reveal the regime of tight coupling. Thanks to the analytic tractability of the discrete-state model, 
the consequences of the transition between the tight- and loose- coupling regimes can be understood for the kinetic 
and thermodynamic properties of the Fi motor. 

An important difference between the continuous and discrete models is that mechanical properties such as the 
friction of the object attached to the 7-shaft in the liquid surrounding the motor or the external torque acting on the 
shaft explicitly appear in the Fokker-Planck equations defining the continuous model, but does not appear in the master 
equation of the discrete model. Indeed, this latter is defined by its transition rates and its rate constants although the 
former also contains a biased diffusive part involving both the friction coefficient £ and the external torque r. In order 
to describe these mechanical properties in the discrete-state model, we thus have to give a dependence on the friction 
coefficient £ and the external torque r in the rate constants, which are also the the constants of the reactions of the Fi 
motor with ATP, ADP, and P;. Accordingly, the discrete-state model is essentially a model of the chemical kinetics 
of the Fi motor with reaction constants fc p (C,T) depending on mechanical parameters. This dependence is fitted to 
numerical simulations of the continuous-angle model. This fitting is performed in the tight-coupling regime where 
the discrete-state model is supposed to correspond to the continuous-angle model. In this way, the two-state model is 
completed by including the dependences on both the chemical and mechanical control parameters. The random time 
evolution of the discrete two-state model can be simulated thanks to Gillespie's numerical algorithm (25l.[26j. 

In the discrete-state model, the tight coupling between mechanics and chemistry implies that the motor is driven out 
of equilibrium by the unique chemomechanical affinity combing the external torque with the chemical potential 
difference of ATP hydrolysis. Thermodynamic considerations shows that the chemical equilibrium is displaced by the 
external torque acting on the motor if this latter strictly functions in the tight-coupling regime. 

Thanks to the solvability of the two-state model, the stationary solutions of the master equation can be exactly 
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deduced, allowing us to show analytically that the mean rotation rate obeys a Michaelis-Menten kinetics with respect 
to the ATP concentration. Moreover, the analytical formula (|38[) is obtained for the rotation rate in the presence of 
both ATP and the products of ATP hydrolysis, i.e., ADP and P ; . 

The highly nonlinear dependence of the mean rotation rate of the 7-shaft (|64[) on the chemomechanical affinity 
(|23p shows that the Fi motor is not functioning in the linear-response regime defined by Onsager's linear-response 
cocfficents, but instead typically runs in a nonlinear-response regime which is more the feature of far-from-equilibrium 
systems than of close-to-equilibrium systems. This remarkable property allows the motor to rotate under physiological 
conditions at about 130 rev/s although it would be a million times slower if the motor was functioning in a linear- 
response regime. 

Furthermore, the crossover between the reaction-limited and friction-limited regimes is well described by the two- 
state model as the friction coefficient £ is increased. Although, the mean rotation rate is nealry independent of the 
friction coefficient in the reaction-limited regime at low friction, it decreases as the inverse of the friction coefficient 
at high friction because the reaction rates (|58|) fitted to the continuous- angle model appropriately take this feature 
into account. 

We have also investigated the behavior of the Fi motor in a surrounding filled with ADP and inorganic phosphate. 
As shown by Eq. (f3"8"]) . the mean rotation rate decreases as the concentrations of ADP and Pi exceed a crossover 
value. The reason is that the release of the products of ATP hydrolysis in the motor is counteracted by the reverse 
reaction of binding of the products to the catalytic sites, which has the effect of slowing down the motor. 

The dependence of the rotation and ATP consumption rates on the external torque is also of great interest because 
it reveals the interval of values of the external torque where the tight-coupling condition holds. It is in this regime that 
the discrete-state model provides a good description of the motor. This comparison shows that the stalling torque 
where the rotation and ATP consumption rates vanish is given in terms of the chemical potential difference of ATP 
hydrolysis and thus depend on the concentrations of the involved species according to Eq. (|67p . 

Moreover, the nonequilibrium thermodynamics of the Fi motor has been developed, allowing us to study the 
chemical and mechanical efficiencies defined in Ref. 27]. In the tight-coupling regime, these efficiencies can reach 
their maximal unit value near the stalling torque. The coupling between mechanics and chemistry is so tight that 
the expected deviations between the stalling torque and the values of the torque where the efficiencies reach their 
maximum value turn out to be very small. The coincidence of these different values of the external torque is the 
consequence of the high degree of tight coupling achieved in the Fi motor, but would not necessarily hold for a motor 
functioning in a loose-coupling regime. 

In conclusion, the two-state model and its comparison with the continuous-angle model of our previous paper fl2| 
is providing a powerful method to study the kinetic and thermodynamic properties of the Fi motor and, especially, 
the coupling between its mechanics and chemistry. 
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APPENDIX A: THERMODYNAMIC RELATIONS BETWEEN AFFINITIES AND CURRENTS 

In this appendix, the nonequilibrium thermodynamics of the Fi motor is presented in the linear regime very close 
to the equilibrium. 



1. The general case 



The molecular motor can be driven out of equilibrium by the two independent affinities which are the mechanical 
and the chemical affinities proportional respectively to the torque r and the chemical potential difference A/i. These 
affinities are the corresponding fluxes or currents can be defined as 

2tt t 

A m = —— «- J m =3V (Al) 



3 k B T 



A c = <-> J c = R (A2) 



19 



in which case the thermodynamic entropy production (|20[) takes the following form: 

-L^ = A m J m + A c J c >0 (A3) 
Kb at 

In general, the currents are nonlinear functions of the affinities 

J m = J m (A m ,A c ) (A4) 
J c = J c (A m) A c ) 

which vanish at the thermodynamic equilibrium where the affinities vanish, A m = A c = 0. Close to equilibrium, the 
currents can be expanded in powers of the affinities, which defines the linear response coefficients as 

J m = L n A m + L 12 A C + 0(2) (A6) 
J c = L 21 A m + L 22 A C + 0{2) (A7) 

The microreversibility implies the Onsager reciprocity relation: 

Lia = L 21 (A8) 

and the non-negativity of the entropy production (|A3D the inequality 

L n L 22 >L 2 12 with L n ,L 22 >0 (A9) 

The coupling between the mechanics and the chemistry is here possible because the motor is attached to a solid 
support and keeps a fixed orientation. The chemistry thus becomes vectorial as well as mechanical. Therefore, there 
is no contradiction with the Curie symmetry principle according to which scalar and vectorial processes cannot be 
coupled together [18l 1 1 Qf| - In isotropic media such as liquids, chemistry is scalar and cannot be coupled to mechanics 
which is vectorial. This coupling becomes possible if the medium is anisotropic as in liquid crystals or in proteins 
attached to the surface of a solid. 

In the plane of the affinities (^4 m ,v4 c ), the curve where the velocity vanishes behaves around the origin as 

J m = W = 0: A c = -^A m + O(Al) (A10) 

while the curves where the ATP consumption rate vanishes is given by 

J c = R = 0: A c = -^A m + 0(A 2 a ) (All) 

-L 22 

Near the origin, the curve (|A10[) has therefore a more negative slope than the curve (|A11|) as the consequence of the 
inequality (|A9[) resulting from the second law of thermodynamics [271 ] . 



2. The case of tight coupling 



In the case of a tight coupling between the mechanics and the chemistry of the molecular motor, these two curves 
have the same slope around the thermodynamic equilibrium point. Indeed, the tight-coupling condition (|2"Tj) reads 

tight coupling: J m = J c = J (A12) 

so that the entropy production becomes 

1 d' S 

tight coupling: l — = AJ > (A13) 

fee dt 

with the chemomechanical affinity ([23|) . which here reads 

A = A m + A c (A14) 

Therefore, the mechanical and chemical affinities are no longer independent in the tight-coupling regime where the 
chemomechanical affinity (|A14[) becomes the unique nonequilibrium driving force. Accordingly, the unique current 
(|A12[) is a function of this unique affinity: 



J = J(A) =LA + 0(A 2 ) 



(A15) 



and all the linear response coefficients become related to each other by 

L = in = L12 = L21 = L22 
In this case, the inequality (1A9|) reaches the equality: 

L11L22 = L\2 

which is also characteristic of tight coupling. 
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(A16) 
(A17) 



3. The efficiencies 

The chemical and mechanical efficiencies (|68l) - (joT))) can be expressed as follows in terms of the affinities and the 
currents: 



Ac Jr. 



1 



Vc = 



A m J m Tjxa 



(Ai8) 



In the regime of linear response where the currents are the linear functions (|A6|) - (|A7|) of the affinities, the efficiencies 
can be written in the form 



Vc = 



1 — e a + 1 



1 



a a + 1 — e rj n 



in terms of the coefficient 



and the constant 



a 



L12 A m 
L22 A c 



e = 1 



T 2 
L11L22 



(A19) 



(A20) 



(A21) 



such that < e < 1 by the inequality (|A9|) . Equations (|A19[) and (| A20|) show that the efficiencies only depend on the 
ratio of the affinities [2jJ • In their respective domains of variation, the chemical and mechanical efficiencies can reach 
the maximum value 



0<7? D 



1 



1 



< 1 



(A22) 



This happens along different curves in the plane (A m , A c ) of the affinities with respect to the curves where the rotation 
and ATP consumption rates vanish since [2JJ 



f?m=?7max if « = -1 + y/s 

J m = 3V = Q if a = -l + e 

J C = R = Q if a = -l 

Vc = f?max if a = -\-sfi 



(A23) 
(A24) 
(A25) 
(A26) 



In the tight-coupling limit, the coefficient (|A21[) vanishes because of Eq. (|A17|) . e — 0, so that the chemical and 
mechanical efficiencies take the value 



1 



Vc 



1 

a- 



An 



(A27) 



according to Eqs. (|A19|) and (1A16|) . In this case, the four conditions (IA23j) - (IA26j) coincide in a = —1 where the 
chemical and mechanical efficiencies reach their maximal value which is equal to unity, 77 max = 1- Therefore, tight 
coupling favors the optimization of the efficiencies. 
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